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ABSTRACT 


During the summer of 1995, an intensive, joint field study called Shallow 
Water Acoustics in a Random Medium (SWARM '95) was conducted by the Naval 
Research Laboratories (NRL), Woods Hole Oceanographic Institution (WHOD, 
University of Delaware (UD), Applied Physics Laboratory of Johns Hopkins 
University (APL/JHU) and Naval Postgraduate School (NPS), among others, in the 
Mid-Atlantic Bight continental shelf region off the coast of New Jersey. 
Environmental and acoustic sensors were deployed as part of SWARM '95 to 
measure and characterize the non-linear internal waves and their impact on the 
Spatial and temporal coherence of the acoustic transmissions. As part of the 
environmental monitoring network, two bottom-moored, upward looking Acoustic 
Doppler Current Profilers (ADCPs) were deployed. A modal, time-series analysis 
of the data captured by the two ADCPs was performed. Highlights of the results 
reveal that: the generation mechanism, in this case, 1s consistent with the lee-wave 
hypothesis of generation; the phase speed is in good agreement with predicted 
phase speeds of the first baroclinic mode; and, the displacement power spectral 


density is significantly modified when soliton wavepackets are present. 
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I. INTRODUCTION 


A. MOTIVATION AND THESIS OBJECTIVES 

In conducting naval operations near any coast, an 
intimate knowledge of the barotropic tide has long been 
regarded as a minimum requirement to enhance safety. 
Through centuries of study, barotropic tides have become 
sufficiently well understood to allow accurate prediction of 
sea-level heights and tidal currents within littoral 
regions. Unfortunately, the ability ©€0 pmedict tical 
motions throughout the water column above the continental 
slope and shelf is not as well understood. When the 
barotropic tide interacts with complex bathymetry, enhanced 
meteorological forcing and the often-stratified vertical 
density structure of the water mass above the shelf, it can 
couple energy into baroclinic modes. As a result of this 
coupling, internal waves, both linear and non-linear, abound 
in the water above the continental shelf. Of particular 
interest are non-linear internal waves, some of which have 
been described as solitary waves (solitons). Although they 
are internal waves and their predominant impact is as 
disturbances within the water column, solitons may exhibit a 
surface manifestation of intermittent choppy and glassy 
seas. In the water column, the internal disturbances caused 
by non-linear internal waves have been observed to be 


Significant enough to potentially impact naval operations 


above the continental shelf. Consequently, as naval 
operations continue to move into littoral regions, gaining a 
better understanding of non-linear internal waves may be as 
important to safe naval operations as has been determining 
the phase of the barotropic tide. 

The focus of this thesis is on a physical description 
of the non-linear internal wave packets, formed near the 
shelfbreak and traveling shoreward on the Mid-Atlantic Bight 
continental shelf. Specifically, Acoustic Doppler Current 
Profiler (ADCP) measurements, obtained from two moorings 
above the continental shelf, are analyzed in an effort to 
answer the following questions: 

1. What is the generation mechanism? 

2. What are the propagation characteristics? 

3. What are the spectral characteristics? 

The ADCP measurements were conducted as part of the 1995 
Shallow Water Acoustics in a Random Medium experiment 

(SWARM '95). The SWARM '95 experimental region is shown in 
Figure 1, along with the bathymetry and locations of the two 
ADCPs. 

The organization of this thesis is as follows: A 
background on SWARM '95 and a brief review of previous 
observations of non-linear internal waves are presented in 
the remainder of this Chapter. The geometry of the ADCP 
measurements and the data processing method used in this 


study are discussed in Chapter II. Analyses on the 


SWARM Experimental Site 


Latitude (degrees) 
Ww a 
60 $ © 
oO =) on 


GO 
de) 


38.5 





38 
-75 -74.5 -74 -73.5 -73 -72.5 -72 -71.5 -71 
Longitude (degrees) 
Figure 1. Locations of the ADCPs in SWARM '95. 


generation mechanism, propagation characteristics and 
spectral characteristics are presented in Chapter III. 
Conclusions are summarized in Chapter IV. 
B. SWARM ‘95 

Under the sponsorship of the Office Of Naval Research 
(ONR), SWARM '95 took place during 20 July to 12 August 
1995 in the Mid-Atlantic Bight continental shelf region off 
the coast of New Jersey (see Figure 1). The primary 
objective of this experiment was to quantify the interaction 
of 10-2000 Hz acoustic signals with shallow water linear and 


non-linear internal waves. In order to accomplish this 


objective, very high quality physical oceanographic data and 
bottom geoacoustic data were collected along with acoustic 
transmission data. In addition to supporting the acoustic 
measurements, the oceanographic data were of sufficient 
quality to support studies of generation, propagation and 
dissipation of the internal waves within the shelf region. 

SWARM '95 was a large, multi-institution, 
multi-investigator project. Investigators from the Naval 
Research Laboratories (NRL), Woods Hole Oceanographic 
Institution (WHOL), University of Delaware, and the Applied 
Physics Laboratory of Johns Hopkins University (APL/JHU) 
participated directly in the cruise. Investigators from 
Naval Postgraduate School (NPS), University of Rhode Island 
(URI), Scripps Institute of Oceanography (SIO), University 
of Miami, and Northeastern University contributed equipment 
and/or participated in the analysis. NPS contributed two 
Acoustic Doppler Current Profilers (ADCPs) and subsequently 
analyzed the measured time-series. 

The collection of the acoustic and oceanographic data 
during SWARM '95 involved three research vessels working at 
the experimental site, two of which occupied the site for 
the full duration of SWARM. The collection included 
successful deployment and recovery of twenty-five 
oceanographic and/or acoustic moorings. A representation of 
the environmental data collection effort is provided in 


Figure 2. The environmental measurements included those 
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from: two bottom-moored ADCPs, thermistor chains, expendable 
bathythermographs (XBTs), a towed conductivity temperature 
depth (CTD) sensor, CTD casts, a vessel-mounted ADCP along 
with imagery from spaceborne and airborne Synthetic Aperture 
Radar (SAR) systems and shipborne radar imagery. 

C. A REVIEW OF SOLITON OBSERVATIONS 

The name of “SWARM” includes the words “random medium.” 
This refers to the environmental fluctuations caused by 
internal waves, particularly the non-linear internal wave 
packets, present in the continental shelf. These waves have 
been described as solitary waves, or solitons, and will be 
discussed herein. 

Scientific literature is rich in articles describing 
solitons in numerous contexts, including: nuclear physics, 
electron-positron plasmas, optical waveguides and 
communications, charged particle energization, conducting 
polymers and others, including oceanography. Of particular 
interest to this thesis are solitons in shallow water. 

According to Osborne and Burch (1980), surface solitons 
were first observed by John Scott Russell in the 19th 
century and Korteweg-de Vries (KdV) published a theoretical 
article about solitons in 1895. Theoretical papers about 
internal solitons were written in the 1950s and 1960s, 
including those by Long (1953, 1956) and Davis and Acrivos 
(1967). Thereafter, numerous empirical results and 


theoretical considerations have been published. 


Some relevant features regarding internal solitons are: 

1) These waves are usually observed in summer when 
they are trapped in a strong and shallow seasonal 
thermocline. [Farmer and Smith (1980); Loder et al. (1992)] 

2) The solitary waveforms are depressions in the 
mvenocline and OLten havera "Sech-” profile. @nhemesin 
(1983)] When an initial solitary wave propagates, it often 
evolves more rank-ordered solitons and exhibits clear non- 
linear and dispersive features such as higher-than-linear 
group velocity, and a decrease in wavelength and amplitude 
toward the rear of the packet. [Farmer and Smith (1980) ] 

3) Solitary waves propagate shoreward, and are often 
generated by a tidally driven flow over sills, continental 
shelf edges, or other major variation in the underwater 
topography. [Chereskin (1983), Loder et al. (1992)] Internal 
waves have also been observed propagating seaward, but will 
not be discussed here. 

4) Because of the shallowness of the summer seasonal 
thermocline and the large amplitude of the coastal internal 
waves, strong surface expressions of solitons have been 
observed with a variety of remote sensors, including 
photographs and synthetic aperture radar (SAR) from 
satellites and space shuttles. [Apel (1978, 1995), Chereskin 
(dlsytes oi ) 

Figure 3 provides a schematic plan view and profile of 


a soliton packet, with typical scales listed in Table 1. A 


spaceborne synthetic aperture radar (SAR) image of the 
surface expression of actual solitons in the vicinity of the 
Mid-Atlantic Bight continental shelf is presented in 

Figure 4. The solitons are the quasi-linear features 
throughout Figure 4, most strikingly is the band from the 
bottom center to the upper-right corner of the image. 

Naturally occurring, large-amplitude internal solitons 
have been reported in many coastal zones of the world such 
as: the Massachusetts Bay, the New York Bight, Gulf of 
California, Andaman Sea offshore of Thailand, the Australian 
North West Shelf, the Sulu Sea between the Philippines and 
Borneo, off the coast of Portugal, off Hainan Island in the 
South China Sea, off the Strait of Gibraltar in the Alboran 
Sea, the Scotian Shelf off Nova Scotia, the Celtic Sea and 
so on. They have also been observed in lakes. [Zhou et al. 
(1991)] Apel et al. (1975) have also located surface 
manifestations of solitons off southwest Africa and off 
equatorial east Africa in the Indian Ocean, while Munro 
(1982) reports on solitons in the Davis Strait. 

Internal solitons can have significant amplitudes. 
Loder et al. (1992) noted that during towed Conductivity 
Temperature Depth (CTD) observations above the Northern 
Georges Banks, solitons with vertical amplitudes of 10-40 
meters and horizontal scales ranging from several hundred 
meters to several kilometers were found. Loder et al. 


(1992) went on to note that these observations were 


Horizontal view 


<< 5 
gu 
aS 


Surface 





h,, Py 


Y Vertical view 
a hoy P2 





Figure 3 Schematic plan view and profile of soliton waves (from Apel 1995). 


Table 1. Typical Scales for Continental Shelf Solitons (after Apel, 1995) 


L (km) 


Ne (m) 
50-500 
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consistent with numerous observations by other researchers, 
including: Farmer and Smith (1980); Chereskin (1983); 
Sandstrom and Elliot (1984); Holloway (1987), as well as 
others. Holloway (1987) describes a bore (soliton) that 
propagates shoreward which produces up to 60 m drops in the 
isotherms, although the 30 minute sampling interval he 
employed may not accurately depict the velocities and 
structure of these solitons. 

No direct effect on naval operations has been 
documented, although these waves have been cited as 
potentially overwhelming the station-keeping capabilities of 
the dynamic positioning systems of offshore drilling ships 
[Osborne et al. (1978), Munro (1982)]. Consequently, 
solitons may also impact naval operations in the following 
ways: 

1. disruption of moored mine fields, 

2. uncontrolled motions of both enemy mines and 
friendly mine clearing equipment during Mine-Counter- 
Measures operations, 

3. sudden "depth excursions" by special warfare divers 
while approaching/leaving coastal waters while submerged, 

4. sudden "depth excursions" by submerged submarines 
while operating above or near the continental shelf, 

5. disruption of acoustic paths while localizing 
and/or prosecuting enemy submarines in shallow water. 


Considering their pervasiveness, significant amplitude 


11 


and possibly significant tactical effect, some explanation 
of the generation and propagation of these waves is in 


order. 
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II. DATA AND PROCESSING METHODS 


A. ADCP MEASUREMENT GEOMETRY 

As part of the SWARM environmental monitoring network, 
two RD Instruments (RDI) Self-Contained Acoustic Doppler 
Current Profilers (SC-ADCPs) were moored 4 m above the 
seafloor, facing upward, at two separate locations, 16.5 km 
apart, along the experimental track. While one of them was 
placed in 103 m of water at 39° 10.000'N, 72° 47.416'W near 
the shelf-break (hereafter referred to as the "Deep Site"), 
the other one was in 75 m of water at 39° 25.33369'N, 
72° 56.5942'W ("Shallow Site" hereafter). Each ADCP was 
configured to transmit 225 pings per ensemble, with an 
ensemble sampling interval of 90 seconds. The acoustic 
frequency for the pings was 307.2 kHz. The four upward 
projected beams were oriented thirty degrees from the 
vertical and the data were recorded in earth coordinates 
(north, east and up each being positive). Each ADCP hada 4 
m sampling blank directly above the transmitter head, then 
it sensed the water column in 4m long depth bins. With 
these settings of profiling parameters, the current velocity 
measurements were accurate to approximately 1 cm/s. 

Measurements at the Shallow Site were collected in 16 
depth bins centered from 67 m to 7 m below the surface. At 
the Deep Site, the data was obtained in 22 depth bins 


centered from 95 m to 11m below the sea surface. Thus, the 
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resultant ADCP data set comprises time series of depth 
profiles of three-dimensional current velocities with a 
sampling time interval of 90 s and a depth interval of 4 m. 
While the Deep-Site time series encompasses data from Julian 
date (JD) 204 to JD 221, the Shallow-Site time series is 4 
days shorter, containing data between JD 209 and JD 222. 
B. DATA PROCESSING 

The processing of the ADCP data entailed three 
Successive steps. The first step was axis re-orientation 
followed by vertical empirical-orthogonal-mode (EOF) 
decompositions. The final step involved temporal filtering 
of the modal coefficients/amplitudes. The filtering enabled 
the decomposition of the ocean variability of three 
different time scales; synoptic, tidal and buoyancy. The 
understanding of soliton production and propagation may be 
gained by studying the relationships between meteorological 
events, the barotropic and internal tides and buoyancy 
motion. The three successive data processing steps, along 
with some important results, are discussed in the following. 

1. Re-orientation 

For the ease of the oceanographic analysis, a regular 
coordinate system (x,y,Z) with the y axis pointing shoreward 
(1.e., upslope) and matching the line that joins the two 
ADCPs was introduced. With respect to the earth-oriented 
coordinates, the y-direction is 306° True. The x-axis has 


an along-slope orientation pointing to 36° True. The z-axis 


14 


1S positive upward with z=0 being the sea surface. The 
direction re-orientation was completed by projecting the 
measured currents, which were in earth-oriented coordinates, 
onto this (x,y,z) coordinate system. The current vector in 
(x,y,z) 1s denoted hereafter by (u,v,w). 

2. Vertical EOF Decomposition 

The time-series of the three particle velocity 
components, upslope (v), along-slope (u) and vertical (w), 
profiled at the Deep and Shallow Sites, were then subjected 
to vertical EOF decompositions. For each velocity component 
at each site, the data covariance matrix was formed and its 
eigenvalues and eigenvectors computed. Mathematically, the 


decomposition can be cast as a linear-algebra problem: 


(DD ‘Ye, - he (1) 


where Dis the data matrix with each row containing the time 
samples at each depth bin and thus the columns correspond to 
depth profiles, and DD’ is the data covariance matrix with 
D’ denoting the matrix transpose of D. The column 
eigenvectors (e;) of DD’ are the vertical EOFs and the 
eigenvalues (A,’*) denote the relative distribution of the 
kinetic energy in each of the empirical modes of the 
velocity component being processed. After the EOFs are 
computed, the time series of the corresponding modal 


amplitude, m;, can be easily obtained by projecting the data 


onto the mode: 
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m=dD. (2) 


Alternatively, one may choose to project the data onto the 
dynamical (1.e., theoretical) modes. However, it is rather 
appealing to calculate vertical modes independently of 
theory, so that the purely empirical results can be used to 
validate the theory. In the bottom half of Figure 5, all 
the "Significant" Shallow-Site EOF modes are displayed, 
except for the zeroth (1.e., barotropic) mode of horizontal 
current which is approximately uniform. By "Significant", 
it is meant that these modes have root-mean-square 
velocities twice as large as the ADCP measurement error of 1 
cm/s. The bottom left panel shows the lst and 2nd 
baroclinic EOFs of the horizontal current while the bottom 
right panel shows the only significant EOF for displacement 
(w). It is important to point out that the EOFs for u ancy 
are almost identical. Modes observed at the Deep Site are 
Similar to those of the Shallow Site. The eigenvalues also 
indicate that the first three EOFs for u and v encompassed 
over 93% of the total horizontal kinetic energy at both the 
Shallow and Deep Sites. 

For comparison, the predicted vertical structure based 
on theory is displayed on the top half of Figure 5. These 
dynamical modes were calculated using a buoyancy frequency 
profile (see Figure 6) derived from a CTD cast at a location 


between the two ADCPs. The EOFs and dynamical modes for the 
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Dynamical horizontal current modes Dynamical displacement modes 





-2 0 2 





—2 0 2 -2 —1 0 1 2 
Figure 5. Vertical structure of currents (left panels) and displacements (right panels). Dynamical modes 
calculated based on an observed buoyancy profile are shown on the top panels whereas the vertical 
empirical orthogonal functions derived from the ADCP data are shown on the bottom panels for comparison. 


First baroclinic mode is denoted by dashed lines and second mode by dot-dashed lines. The zeroth 
barotropic mode is not shown. 


horizontal currents are strikingly similar. However, as for 
the vertical current, dissimilarities between the EOF and 
the lst baroclinic dynamical mode can be easily observed. 
An explanation is that the EOF for w may actually be made up 
of a few dynamical modes with highly energetic, phase-locked 
soliton wave packets. In support of this explanation is the 


fact that the lst and 2nd horizontal-current EOFs are highly 
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Figure 6. A Brunt-Vaisala frequency profile (left) and a temperature gradient profile (right). Both profiles 
are derived from a CTD cast taken between the two ADCP sites with some smoothing. 


correlated with correlation coefficients larger than 0.7 at 
both the Shallow and Deep Sites. The empirical 
decomposition was able to discern these modes in u and v 
because of the more energetic tidal motions. But the tides, 
including the internal ones, have little signature inw. It 
was also found that an almost perfect fit to the EOF for w 
can be obtained using only the first three dynamical modes, 
with the first baroclinic mode being the dominant mode in 
SuUCh camEIts 

3. Temporal Filtering 

The velocity modal-amplitude time-series were then 
filtered using lowpass, bandpass and highpass Butterworth 
filters, thereby dissecting the motions into three distinct 
frequency bands. These three frequency bands are hereafter 


referred to ask 
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1. "synoptic," with periods longer thane teahourc, 

2. “tidal," with periods between 5 and 18 hours, and 

3. "buoyancy," with periods shorter than 5 hours. 

The synoptic band should include oscillations directly 
related to the synoptic weather pattern as well as inertial 
waves. The tidal band should consist of semi-diurnal tides. 
The internal waves including the non-linear solitons are 
expected to reside in the buoyancy band. Figures 7 and 8 
show some of the filtered Deep-Site modal horizontal 
currents in the tidal and buoyancy bands, for two recording 
meriods, JD 208 through JD 21235 and JD) 212.5 theeugh 
JD 217, respectively. The observed relations between the 
large-amplitude solitary wave packets, the barotropic and 
internal tides, and the synoptic atmospheric forcing are to 


fPemcaiscusseaq in thegnext secti#en. 
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Figure 7. A segment of the observed Deep-Site modal amplitudes between JD 208 and 212.5, showing (top) the upslope tidal currents of the barotropic 
(solid line) and the 1st (dashed line) and 2nd internal (dot-dashed) modes, (middle) the upslope buoyancy currents of the Ist mode, and (bottom) the along- 
slope buoyancy currents of the 1st mode, respectively. In the middle panel, three of the soliton packets are labeled as A, B and C. 
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IIT. ANALYSIS 


A. GENERATION 

The ADCP data reveal that there were multiple soliton 
generation sites in the vicinity of the experimental track, 
as one might surmise after viewing satellite SAR imagery, 
Figure 4. From the u and v buoyancy-band modal time-series, 
as displayed in the bottom two panels of Figures 7 and 8, it 
will be shown that one of the dominant sites of generation 
was directly off the shelfbreak adjacent to the Deep-Site 
ADCP. The solitons generated near the Deep-Site ADCP 
propagated directly upslope (i.e., in the y direction) 
showing significant signatures in the v component but not in 
the u component. Examples of these locally generated 
soliton packets are denoted by "A" and "C" in Figure 7. In 
the same figure, an example soliton packet that was not 
generated near the Deep-Site ADCP, but was generated to the 
southwest of the experimental region, is labeled "B". This 
"B" packet had a u component of comparable speed to that of 
the v component. The overall environmental condition was 
extremely complex, with distantly generated soliton packets 
also influencing the local condition, either as isolated 
events or by terrerine with the locally generated packets. 
This complexity is more apparent in the Shallow-Site ADCP 
data as the distantly generated solitons tended to have 


larger amplitudes there. Those distantly generated soliton 
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packets appeared to come primarily from the southwest and 
the south-southwest. Further analysis in conjunction with 
the SAR data could help to localize where these distant 
Sites are. 

The top panel of Figure 7 depicts a period when the 
barotropic and internal tides are in-phase in the upslope 
direction at the Deep-Site, at least until approximately 
JD 211.5. This was a period when the distantly generated 
solitons were weak and sometimes absent and the locally 
generated solitons predominated. These soliton packets rode 
on the leading edge of the internal tides and coincide with 
the onset of the barotropic tide. However, as shown in 
Figures 7 and 8, between JD 211.5 and JD 214.5, the 
distantly generated solitons appeared more frequently and 
caused considerable interferences. The distantly generated 
solitons became difficult to isolate from the locally 
generated ones. It was also observed that, in this period, 
the barotropic and internal tides are not exactly in phase. 
This 1s probably due to the interferences between the 
internal tides, upon which the solitons, emanating from the 
multiple source sites, initially rode. 

The generation of solitons in the region of the shelf- 
break may be explained by different, although not 
necessarily competing, hypotheses and observations. Apel et 
al. (1995) suggested two possible mechanisms for the 


generation: coupling of energy between tidal modes and the 
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release of a standing lee wave off the shelfbreak. In the 
former generation mechanism, energy from the barotropic tide 
is scattered into the internal modes within the tidal band 
at the shelf break as the mode angle of the internal 
propagation characteristics coincides with the angle of the 
shelf slope. The internal tides then scatter energy non- 
linearly into the buoyancy band as they propagate shoreward. 
For reasons explained below, this is not expected to be the 
generation mechanism. 

The latter mechanism may explain observations of 
soliton-like internal waves by Loder et al. (1992) and 
Brickman and Loder (1993) in the region of the Northern 
Georges Bank shelf break. In this hypothesis, a steady 
offshore flow, usually provided by the ebbing barotropic 
tide, is required to establish a spatially oscillating 
depression (i.e., lee waves) in the pycnocline directly 
offshore from the shelf break. These depresSions are a 
result of the sudden increase of the relative VOLLIG Hey 
caused by the stretching of the water column. The lee wave 
1s eventually released from its phase locked position when 
the flow reverses. At the onset of the flow reversal, which 
typically occurs at slack tide, the lee waves, then, travel 
up the shelf as internal tides coupling energy into the 
buoyancy band as soliton wave packets. In the absence of 
Ekman flow, both generation mechanisms would produce common 


Signatures in the ADCP observation, in that soliton packets 


Piss 


would be generated at the onset of the flood cycle of the 
barotropic tide, and that they would be riding on the 
leading edge of the internal tides propagating upslope. 
However, as described below, a synoptic wind event existed 
during the experiment, allowing us to discard the former 
hypothesis. 

In Figure 9, the upslope component of the synoptic 
barotropic current 1S superimposed onto the upslope tidal 
current. The synoptic current shows that a weather pattern 
started to move in or develop at approximately JD 214, and 
is confirmed by wind speed readings provided in Figure 10. 
This storm event introduced a synoptic shoreward flow which 
gradually ramped to a peak in excess of 0.3 m/s. This 
synoptic event suppressed and eventually negated any 
offshore flow through the end of the ADCP recording. Of 
note is that during this late period, no solitons were 
locally generated despite strong tidal flow. The 
implication is that the generation mechanism is consistent 
with the lee-wave hypothesis. 

B. DYNAMICAL MODAL RELATIONS 

Central to the subsequent analyses of the ADCP data on 
soliton propagation and spectral characteristics is the 
establishment of the relations between modal displacements 
and modal horizontal and vertical particle velocities. 
These important relations, along with some theoretical 


vertical-structure calculations, are summarized in this 
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Figure 9. Upslope synoptic currents (dashed line) with tidal currents superimposed (solid line), showing 
the onset of a synoptic storm event at approximately JD 214 that prevents any offshore flows in the later 
period. The time-series shown here correspond to the modal amplitudes of the barotropic mode. 
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section. 
The particle velocities (u, v, w) and displacements 
(Nn) associated with the ocean variability can be decomposed 


as linear combinations of vertical modes: 


(uv)=)/ (u (ay), v yt] @ F @) G2 
w=))w (xy,t) W (2) (3b) 
n=do 1 Gay.) W,@) (3c) 


where F, and W, are the vertical modes of horizontal motion 
and vertical displacement or motion, respectively, q, arethe 
normalization constants for F,, and u,, v,, w, and n, are the 
corresponding modal coefficients (i.e., amplitudes) 
describing the temporal and horizontal structure. For 


internal disturbances, the vertical-displacement modes W, 
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Figure 10. Meteorological data for Julian Date 221 showing wind speeds from land 
stations and buoys in the vicinity of the experimental site. The buoy reports display 
winds of 20 knots generally from the NNE and NE. 


are governed by the Sturm-Liouville equation: 


2 
N (2) 
de? oe e as 0 (4a) 


n 





subject to the rigid-lid boundary conditions: 


W (0) = W(-H) = 0 (4b) 
where N is the buoyancy frequency, c, are the eigenvalues 
which are the modal phase speeds of linear waves in a 


nonrotating fluid, i.e., for wave frequency much higher than 


the inertial frequency, and H is the water depth. 
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The vertical modes of horizontal and vertical motions 


are related by 


dw 


n 


F SS 
n dz 


and according to the Sturm-Liouville conditions, they 


satisfy the following orthogonal conditions: 


0 
i NW W dz = 0; men (6a) 
-H 

and 
0 
[FF dz = 0; men. (6b) 
a n im 


Our choices for normalization are 
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] 2 

a | ¥,@ = ] (7a) 
-H 

and 
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1 2 1 

al F.% a Tae (7b) 
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With such normalizations, the modal coefficients have the 
physical meanings of phase-preserved, root depth-averaged 
square quantities. The 1st and 2nd dynamical baroclinic 


modes of horizontal current (1.e., F, and F,) and of 
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vertical displacement or current (1.e., W, and W,) are shown 
on the top panels of Figure 5. These theoretical vertical 
structures were calculated based on a buoyancy profile 
observed by a hydrographic station between the two ADCP 
Sites. This buoyancy profile along with the associated 
temperature gradient profile are shown in Figure 6. It is 
seen that the buoyancy frequency ranges from approximately 5 
cycles/hr near the bottom to a maximum of 30 cycles/hr in 
the thermocline with a columnar average of about 10 
cycles/hr. It is also seen that the temperature gradient 
can be as large as 0.7 °C/m at the thermocline, implying 
that a 10-m internal wave can cause a 7 °C temperature 
change and a 15-m wave can produce temperature fluctuations 
as much as 10.5 °C. Such large fluctuations are often seen 
in the thermistor-chain records during the passage of 
soliton wavepackets. [Lynch, personal communication] 

Using the modal expansion in the fundamental 
hydrodynamics equations, neglecting rotational effects and 
energy exchange (i.e., coupling) between modes, and assuming 
weakly nonlinear, plane progressive waves propagating ina 
specific direction, Ostroskiy (1978) and Apel et al. (1995) 
have shown that the modal displacements h, are governed by 


the KdV equation: 
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where y' is the direction of propagation of the wavefront, 








and 
3c > 0 
a= fF ae (8b) 
n 2H F n 
and 
‘e q° 
B _ fon (8c) 
n 2 


are environmental parameters accounting for the effects of 
buoyancy, current, current shear and water depth. To 3rd- 
order, the modal displacements and particle velocities obey 
the relations (Ostroskiy, 1978): 

on ' ov Ss a(n wy) 


+ q ——— = 0 (9a) 
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with 


q. 0 


s = [Fd (9b) 
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where v', is the particle velocity parallel to y', and 
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The parameters s, and k, express the degree of non-linearity 
in the relations between n, and v', and n, and w,, 
respectively. 
C. SOLITON PROPAGATION 

For small amplitude modal displacements, the 3rd and 
Ath terms of the KdV equation (9) are negligible. In such 
case, the solution is that of the linear internal waves 
_ propagating at phase speeds of c,. However, for 
displacements having large enough amplitudes and steepness, 


the KdV equation admits the well-known soliton solution 
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where n,, 1S the peak modal amplitude of the soliton. The 
superscript (s) is used hereafter to denote modal amplitudes 
associated with solitons. The soliton has a nonlinear 


characteristic width: 


128. 


A = (1 1b) 
n a T 
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and propagates with a nonlinear phase speed 
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(llc) 


Both the soliton characteristic width and phase speed depend 
on the soliton peak displacement. The characteristic width 
is inversely proportional to the square root of f,,, whereas 
the phase speed is linearly proportional to n,,. The 
implication is that the larger the peak displacement, the 
faster the soliton propagates and the narrower, or steeper, 
the soliton is. 

Using (9) and (10) and the soliton displacement 
solution, one can easily obtain the corresponding soliton 
modal horizontal and vertical velocities. Neglecting 3rd 


and higher order nonlinear terms, they are: 
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Using these soliton modal relations with environmental and 
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non-linear parameter values calculated based on the observed 
buoyancy frequency profile, predictions on the evolution of 
solitons can be attempted. To illustrate the predicted 
temporal structure of a soliton and the dramatic effect of 
non-linearity, a soliton solution with a 10-m peak modal 
displacement along with the corresponding modal horizontal 
and vertical velocities are displayed in Figure 11. Of 
particular interest is the dot-dashed curve in the first 
panel. This curve represents the linearized solution 
associated with the particle velocities displayed in the 
lower panels of Figure 11. Clearly, the non-linear 
advective effect is important as it significantly enhances 
the displacement. The modal soliton solutions for the first 
two baroclinic modes are always positive. The positive sign 
is entirely consistent with the data. In view of the 
corresponding mode functions shown in Figure 5, the positive 
Sign implies that the solitons in Mode 1 cause only 
depressions of the isotherms over the entire water column. 
Whereas for Mode 2, it is a depression below 20 m. 

Figure 12 shows the predicted dependence between soliton 
peak modal horizontal current, peak modal vertical current 
and peak displacement. Again, the predicted relation 
between the soliton peak horizontal and vertical currents is 
rather consistent with the data. The observed peaks of the 
Mode 1 soliton vertical current range from 2 cm/s to 10 cm/s 


whereas the corresponding horizontal currents have 
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Figure 11. The evolution of modal displacement (top), horizontal current (middle) and vertical current 
(bottom) associated with a single soliton. The soliton temporal structure (solid lines) was calculated based 
on weakly non-linear theory. The dot-dashed line in the top panel represents the result of linearization, 
showing large error. 


peaks between 5 cm/s and 30 cm/s. It 1s interesting to 
point out that the caps of these observed current peaks 
coincide roughly to where the predicted curve of peak 
vertical current versus displacement turns (see bottom panel 
of Figure 12). This turning of the curve corresponds to 
where the weakly non-linear soliton theory starts to fail 


and where the third- and higher order terms become 
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Figure 12. Theoretical relations between soliton peak displacement and (top) peak horizontal current, and 
(bottom) peak vertical current for the first (dashed lines) and second (dot-dashed) baroclinic modes. 


important. Near this cap, the solitons are extremely steep 
(l.e., narrow) with displacements exceeding 15 m. It is not 
unlikely that they would break and contribute significantly 
to vertical mixing. 

The alignment of the two ADCPs with the adjacent 
soliton generation site has allowed for a verification of 
the theoretical soliton phase speed as passage times at two 


ADCPs which were separated by a known distance could be 
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observed. Before the onset of the synoptic storm event 
which mitigated the production of solitons, there was a 
Gommon. time period, JD 209.5 through JD 216/57 82m whten corn 
ADCPs were recording simultaneously. Using this five-day 
segment of the time-series, it was possible to track a total 
of eight non-interfered soliton wavepackets in the first EOF 
that were generated off the immediate shelfbreak (1.e., with 
relatively small u velocities). While the travel time from 
the Deep Site to the Shallow Site of the leading soliton in 
each packet provided an estimate of its phase speed, its 
peak vertical velocity provided estimate of its peak 
displacement. The results are marked by circles on the top 
panel of Figure 13, in which the theoretical phase-speed 
curves for the ist and 2nd dynamical modes are also 
displayed. Qualitative agreement is noted between the 
observed phase speeds and the predicted phase speeds of the 
1st baroclinic mode. This is not surprising because the EOF 
for displacement (w) is mostly matched by the list dynamical 
mode. A real surprise, however, is the high correlation 
(above 0.7) between the modal amplitudes of the lst and 2nd 
EOFs of the buoyancy-band upslope horizontal current (v), at 
both the Deep and Shallow Sites. High correlation was only 
anticipated at the Deep Site but not the Shallow Site, which 
was more than 16 km away from where the soliton packets were 
generated. The implication is that Mode 2 did not propagate 


with the much smaller phase speed in accord with the weakly 
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Figure 13. Theoretical relations between soliton phase speed (top), 50% width (middle) and 50% time 
duration (bottom) and peak modal displacement for the first (dashed lines) and second (dot-dashed lines) 


baroclinic modes. Observed soliton phase speeds are indicated by circles showing agreement with the 
predicted first-mode phase speed. 


nonlinear, de-coupled soliton theory. Instead, it 
propagated together with Mode 1. A phase-locked mechanism 
apparently existed. Therefore, an outstanding question is 
what causes this phase-locked behavior. Could it be due to 
mode coupling such that Mode 1 replenished the Mode 2 energy 
along the way? Could it, instead, be due to higher-order 
non-linearity, since the phase speeds and particle speeds of 


Mode 2 have the same magnitude? Or, lastly, might it be due 
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to tidal shear which may have strongly modified the tidal 
band vertical mode structure? These questions remain open 
for future investigations. 

In the bottom two panels of Figure 13, the predicted 
50% width and time duration, for both Modes 1 and 2, as a 
function of the soliton peak modal displacement are plotted. 
We use "50%" to denote where, or when, the displacement 
drops to 50% of its peak value. The theoretical curves for 
Fede i suggested that the observed solitons should have 50% 
widths of about 100 to 200 m and 50% time durations of about 
100 to 300 s. These width and duration estimates are 
consistent with observations. 
D. DISPLACEMENT POWER SPECTRAL DENSITY 

The displacements of isotherms are directly related to 
sound-speed perturbations. Therefore, for underwater 
acousticians, the power spectral density of displacements is 
a more preferable characterization of the acoustic random 
medium than those of the particle velocities. In this 
Spirit, an attempt to provide some estimates and discussion 
of the displacement power spectral density will be presented 
next. 

Displacement time-series can be derived from the 
vertical-current time-series. The relation between modal 
vertical velocity and modal displacement as stated in (13) 


can be integrated to give 


Sh, 


n- km. = fw de. (15) 


It 1s noted that this integral relation for solitons is also 
valid for linear internal waves residing in the energy- 
containing frequency band of the soliton packets. For these 
linear waves, the nonlinear k,n,’ term is much smaller than 
NH, and hence contributes little in the overall equation. 

The soliton energy-containing band can be easily estimated 
by performing a Fast Fourier Transform (FFT) on just a few 
individual soliton packets. Such a band is 40 cycles/day 
and above. Therefore, if it 1s assumed that the greater- 
than-40-cycles/day ocean variability consists of linear 
internal waves and solitons only, (15) can be used to 
convert the observed vertical currents into displacements in 
this band. 

Thus, the procedure in obtaining the modal displacement 
time-series involved, as a first step, filtering the 
observed modal vertical velocity into a lowpass time-series 
and a highpass time-series, both with a cutoff at 40 
cycles/day. After integrating both series in time, the 
highpass one was then corrected, based on (15) through 
polynomial root finding. Finally, the lowpass and corrected 
highpass displacement series were added to give the overall 
displacement time-series. Partial segments of the 
corrected, highpass modal displacements at both the Shallow 


and Deep Sites are shown in Figure 14. Having a common time 
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window, the two segments show the evolution of four soliton 
packets over two tidal cycles. Although it may be difficult 
to see without the aid of the horizontal currents, there 
were actually two packets per cycle with the locally 
generated packet leading. 

With the expectation that the sel characteristics 
would be dissimilar between time periods with little soliton 
activities and those with heavy activities, it is only 
logical to treat the active and inactive periods separately 
in the spectral estimate. The displacement time series were 
carefully cut into many non-overlapping sections, each six 
hours long. The six-hour window was chosen because each 
soliton packet was typically four to six hours long. For 
each ADCP, the six-hour sections were then assembled into 
two groups, one with soliton packets and the other without 
soliton packets. The spectral density of each group was 
then estimated through ensemble averaging of the magnitude 
squared sectional FFTs. 

Figure 15 shows the resultant power spectral density 
estimates of the modal displacements, at the Shallow and 
Deep Sites, for periods with and without soliton activities. 
For the solitonless periods, both the Shallow and Deep-Site 
spectra exhibit a frequency dependence of approximately f°” 
over the band from 10 to 150 cycles/day. But beyond 150 
cycles/day (1.e., the buoyancy frequency), the solitonless 


spectrum falls off much more gradually in the Shallow Site 
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Figure 14. A segment of highpass time-series of modal displacements, at the Shallow (top) and Deep 
(bottom) Sites. The upslope barotropic tidal current (dot-dashed line in units of cm/s) is superimposed to 
show that this segment encompasses two tidal cycles with strong soliton packets generated near the Deep 
Site during the tidal flood. 


than in the Deep Site. The f°?’ dependence below 150 
cycles/day is rather remarkable in that the same dependence 
was observed by Lynch et al. (1996) in the Barents Sea. 
Since the Barents Sea 1s also on a large shelf, but at much 
higher latitude, the commonalityeof this fw" dependence fer 
shelf regions 1s an interesting question to ask. 


It is apparent from Figure 15 that the spectral effects 
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Figure 15. Power spectral density estimates for periods with no soliton activities (solid curves) and periods 
of strong soliton activities (dot-dashed curves) for both the Shallow (left panel) and the Deep (right panel) 
Sites. The solid straight lines represent linear least-squares fit to the solitonless logarithmic spectra in two 
different frequency bands. The slopes of these fits are shown in parentheses. 


of the soliton packets include the elevation and whitening 
of the band from 40 to 200 cycles/day with slight resonance 
toward the higher end of this band. This slight resonance 
1s associated with the full durations of the individual 
solitons. For peak modal displacements of 8 to 10 m, the 
soliton full durations are predicted by theory to be 
approximately between 600 and 500 s. The reciprocals of 


these durations give 140 to 170 cycles/day. The spectral 
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elevation is significant, revealing an order-of-magnitude 
increase of potential energy in this elevated band. It was 
also observed that the soliton-active spectral density falls 
off much more rapidly beyond 200 cycles/day at the Shallow 
Site than at the Deep Site. This difference might very well 
be due to the breaking of the high-amplitude solitons as 


they propagated shoreward. 


IV. CONCLUSIONS 


Analysis of the nearly 3 week-long ADCP records 
obtained during the summer near the shelf break in the Mid- 
Atlantic Bight reveals that non-linear internal waves are 
evident and appear to be generated at multiple source sites 
in the vicinity of the SWARM '95 experimental site. It was 
observed that generation of non-linear internal waves ceased 
when the off-shelf tidal flow was countered by atmospheric 
forcing during the later part of the experiment. The 
atmospheric forcing induced a synoptic shoreward flow which 
negated any offshore flow coincident with a cessation of 
soliton generation near the experimental site. These 
observations are consistent with the lee-wave hypothesis for 
soliton generation. 

A further finding from the time-series was that the 
observed phase speed of the solitons was in good agreement 
with predicted phase speeds of the first baroclinic mode. 
This was not surprising as the EOF for displacement is 
mostly matched by the first dynamical mode. What was 
surprising was the strong correlation between the modal 
amplitudes of the lst and 2nd EOFs of the buoyancy band 
upslope horizontal current (v) at both the Deep and Shallow 
Sites. This implied that the 2nd EOF propagated phase- 
locked with the lst EOF. 


Lastly, the displacement power spectral density 
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estimates of the soliton-active periods and the solitonless 
periods provided a look at the potential energy associated 
with solitons. First, an increase of potential energy in 
the soliton-active spectrum is apparent when compared with 
the solitonless spectrum. Another point is that the 
soliton-active spectral density falls off much more rapidly 
beyond 200 cycles/day at the Shallow-Site than at the Deep- 
Site. This difference may be caused by breaking of the high 
amplitude solitons as they propagate upslope. 

In order to further the understanding and 
characterization of solitons in SWARM '95, it will be 
necessary to include the other, complementary, data sets in 
future analyses. Hydrodynamic models of the wave packets 
may serve as a framework for a joint analysis of the 
different observations. The study of the impact of solitons 
on acoustic signals may then be enhanced. A further goal 
will be to, eventually, predict soliton generation, 
propagation and dissipation. Hopefully, this will lead to 


predicting the impact of solitons on naval operations. 
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